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We add an energy-independent Hamiltonian to the standard flavour oscillation one. This 
kind of physics might appear in theories where neutrinos couple differently to a plausible non- 
zero torsion of the gravitational field or more dramatically in the presence of CPT-violating 
physics in the flavour oscillations. If this contribution exists, experiments at higher energies are 
more sensitive to their free parameters, and flavour conversion could be severely modifled. We 
show that this new physics modifies the neutrino mixing angles and find expressions that relate 
the new, effective, angles to the standard oscillation parameters Am^j, 9ij and Sep and to the 
parameters in the new-physics Hamiltonian, within a three-neutrino formalism. We consider 
scenarios where the new parameters allow for extreme deviations of the expected neutrino flavour 
ratios at Earth from their standard values. We show that large departures of the standard flavour 
scenario are plausible, which would be a strong hint of the violation of a conserved symmetry. 



1. Motivation and problem 

By now it has been demonstrated that neutrinos can change flavour. Experiments with 
solar, atmospheric, reactor and accelerator neutrinos have established that there is a nonzero 
probability that a neutrino created with a certain flavour is detected with a different one after 
having propagated for some distance, and that this probability is dependent on the propagated 
distance, L, and the neutrino energy E. 

The standard mechanism that explains neutrino flavour transitions makes use of two different 
bases: the basis of neutrino mass eigenstates, which have well-defined masses, and the basis of 
neutrino interaction states -the flavour basis- which are the ones that take part in weak processes 
such as W decay. The two bases are connected by means of a unitary transformation U, so that 
we can write each of the flavour states as a linear combination of the mass eigenstates 
i.e. 



where the coefficients Uai are components of the unitary mixing matrix that represents the 
transformation. Assuming the existence of three active neutrino families (a = e,/i, r), as 
indicated by experiment [Ij, the summation index runs i = 1, 2, 3 and C/ is a 3 x 3 matrix. Using 
this definition for the flavour states it is straightforward to conclude that flavour transitions occur 
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because neutrinos are massive, and because different mass eigenstates fiave different masses. Tfiis 
can be seen, within tlie tliree-neutrino formalism, by writing tlie transition probabilities from 

fa to Uf^ as 

= 5^,-4^Re (4) sin^ (^A + 2 J] Im sm (^l) , (2) 

i>j \ / i>j \ / 

where Am?^ = mf — m^, with rrii the mass of the i-th eigenstate, and J^^ = U*^UpiUajUpj. 
Clearly, if Am?^ = 0, no transitions occur. Note the 1/E dependence on the energy associated 
with this standard, mass-driven, oscillation mechanism. 

The experiments that have studied neutrino flavour transitions [2] have been designed to 
detect neutrinos with energies that range from a few MeV (solar neutrinos) to the TeV scale 
(atmospheric neutrinos). The mass-driven mechanism mentioned in the preceding paragraph 
has been experimentally confirmed within this energy range. Notably, data from the Super- 
Kamiokande atmospheric neutrino experiment [3] was used to find an energy dependence of the 
oscillation probability of -E", with n = —0.9 it 0.4 at 90% confidence level, thus confirming the 
dominance of the mass-driven mechanism in this energy range, and relegating any other possible 
mechanism due to new physics to subdominance. It is possible, however, that one or more of 
such subdominant mechanisms become important at higher energies. 

In the present paper, we have explored a possible scenario where there is an extra oscillation 
mechanism present which results in an energy- independent contribution to neutrino oscillations. 
This mechanism, though subdominant in the MeV- TeV range, might become dominant at higher 
energies. We will deferr the detailed treatment of how the contribution is introduced to Section 
[2] and focus now on the possible mechanisms behind it. 

An energy-independent contribution to neutrino oscillations is introduced by a vector coupling 
of the form 

£ = V^bf-f^'i^'^ . (3) 

Such a term results in an energy-independent oscillation phase bij = hi — bj, with the bi 
eigenvalues of the b matrix, so that the phase in Eq. ([2]) becomes 

Am^j Amfj 

and the mixing matrix U is modified as well. Specially relevant to the detection prospects of 
new physics is the fact that not only the oscillation phase -which, as we will see, is lost when 
considering neutrinos from distant astrophysical sources-, but also its amplitude, is modified by 
such an energy-independent contribution. 

A vector coupling like that of Eq. ([3]) could be induced by violation of Lorentz invariance 
which in turn results in a violation of CPT symmetry. A realisation of this possibility is the 
Standard Model Extension [Sj |6], which includes CPT-odd terms in the Lagrangian, such as 
that of Eq. ([3]) . A different possible mechanism is a nonuniversal coupling to a torsion field [4j . 

In the next Section, we will argue that a possible energy-independent contribution to flavour 
oscillations could become important, with respect to the standard oscillation terms, when the 
neutrino energy is high (PeV or more). The highest-energy available fiux of neutrinos is the 
expected ultra-high-energy (UHE) flux from astrophysical sources -notably, active galaxies- 
which are located at distances in the order of tens or hundres of Mpc. For these conditions, the 
flavour-transition probability in Eq. ([2]) oscillates very rapidly, and so we are moved to use its 
average value instead, which is obtained by averaging the oscillating terms in the expression, i.e. 

(i^a/j) = J]|f/aiP|C//3d' • (5) 



Thus the information in the oscillation phase is lost, but the modifications due to the energy- 
independent new physics is preserved in the elements of the mixing matrix. 



2. Theoretical framework 

2.1. Standard mass- driven flavour oscillations 

As we saw in the previous Section, the standard scenario of neutrino oscillations considers each 
neutrino of a given flavour to be a linear combination of neutrino mass eigenstantes, namely 

i 

where we have renamed the mixing matrix Uq. The mass eigenstates satisfy Schrodinger's 
equation and so propagate as 



where we have assumed that <C E, so that p = ^^ E"^ — mf ~ m^/{2E), and that, because 

neutrinos are highly relativistic particles, t ~ L (in natural units). Thus, after having propagated 
for a distance L, the neutrino created with flavour a has become 

K(L)) = j;[C/oLe-i*^|...) . (8) 

i 

We see that each mass eigenstate acquires a phase that depends on the value of its mass: the 
interference of these phases is the source of the flavour transitions. At detection, the probability 
amplitude of seeing a is 

2 2 
« ^ .m _ _ .m- 

id i 

and so the probability Pap = {{I'pli^a for the transition Ua —^ turns out to be Eq. ([2]). 

(The reader is referred to [7J for a more extended treatment of the basics of neutrino oscillations.) 

The oscillation Hamiltonian in the previous description of standard neutrino oscillations was 
written in the basis of mass eigenstates. An alternative, but equivalent, description of the same 
phenomenon can be achieved by writing the Hamiltonian in the flavour basis instead. Using a 
Hamiltonian in the flavour basis will allow us to introduce contributions from new physics in 
a more straightforward manner. If a neutrino is produced with flavour a, then, after having 
propagated for a distance L, its evolved state will be 

\u^{L)) = e-''"-'^\u^) , (10) 

where the oscillation Hamiltonian Hm is the one corresponding to the standard, mass-driven, 
mechanism, and is written in the flavour basis. Hm is related to the Hamiltonian in the mass 
basis -the "mass matrix"- through a similarity transformation that makes use of the unitary 
mixing matrix Uq: 

t diag (O, Am^i , Am?i) . ^ ^ 

Hm = U^HUl = Uo ' '-^U^o . (11) 

Uq is the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix, which can be written in terms of 
three mixing angles, 612, 0i3 and 623, and one CP-violation phase, 5cp, as 

C12C13 S12C13 5136"*'' 

Uo{{Oij} ,Scp) = \ -S12C23 - ci2S23Si3e*'' C12C23 - si2S23Si3e"' S23C13 

S12S23 - Cl2C23Sl3e*'' -C12S23 - 5120235136"^ C23C13 



with Cij = cos {9ij), Sij = sin (Oij). 

Using the latest data from solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K and MINOS) experiments, the authors of [8j found the best-fit values (Icr) 
of the standard oscillation parameters to be 

Am^i = 7.65t[}:^;j X 10"^ eV^ , |Amii| = 2.40t|]:J? x 10"^ eV^ (12) 

sin2 (^12) = 0.304l0:°22 , sin2(%)=0.0ltnn , sin^ (^23) = CSO^^e • (13) 

The are no experimental values for 6c p presently. 

2.2. Adding an energy-independent Hamiltonian 

Motivated by the vector coupling considered in Eq. ([3]) , and in analogy to the standard oscillation 
scenario, we can introduce an energy-independent contribution in the form of the Hamiltonian 
(also in the flavour basis) 

Ht = Ubdmg{0,b2i,bsi)ul , (14) 
where bij = bi — bj. Following [9], we write the mixing matrix in this case as 

U, = diag (0, e^<^^ e^<^^) C/q {{Ouj} , 5,) . (15) 

The mixing angles associated with this Hamiltonian are ^612, ^bis, ^623) and 5h fills the role of 
5c p in the standard Hamiltonian. The two extra phases, (j)2 and (^3, appear because of how the 
mixing matrix between flavour and mass eigenstates was defined in Section [2.11 

H[, is dependent on eight parameters -two eigenvalues (621, ^31)1 three mixing angles (0fei2, 
^613, db2z) and three phases {5b, 4>2, 03 )- whose values are currently unknown. Experimental 
upper limits |9], however, have been set for 621, using solar and Super-Kamiokande data, and 
632, using atmospheric and K2K data: 

621 < 1.6 X 10^2^ GeV , 632 < 5.0 x 10'^^ GeV . (16) 

The full Hamiltonian, including standard oscillations and the energy-independent 
contribution, is then 

Hf = Hm + Hb. (17) 

In Section [H we saw that Hm has been experimentally demonstrated to be the dominant 
contribution to the oscillations in the low energy (MeV-TeV) regime: according to Eq. (I16p . 
the values of the bij are too low for the new physics, if there is any, to manifest at these energies. 
The 1/E dependence of Hm, however, allows us to explore the possibility that, at higher energies, 
when its contribution is reduced, the effect of an energy-independent Hamiltonian Hh becomes 
comparable to it or even dominant. Such energy requirement is expected to be fuUfilled by the 
UHE astrophysical neutrino flux (see Section [TJ) 

We would like to write the flavour transition probability corresponding to this Hamiltonian 
in a form analogous to Eq. @. In order to do this, we need to know what is the mixing matrix 
Uf between the flavour basis and the basis in which Hf is diagonal. According to basic linear 
algebra, this is achieved simply by diagonalising Hf, finding its normalised eigenvectors, and 
building Uj by arranging them in column form. The components of the resulting matrix are 

in general complicated functions of the standard mixing parameters ({%}, |Amfj|, 5cp) and 

of the parameters of Hi, {{Otn.j}, {hj}, 5i,, (j)2, 4'3)- In analogy to Eq. the average flavour 
transition probability associated to the full Hamiltonian Hf is then 

{P^^) = E\^Uflf\[Uf],f ■ (18) 

i 



By comparing the mixing matrix obtained by diagonalisation of Hf with a general PMNS matrix, 
Eq. ()12p with mixing angles @ij and phase Sf, we are then able to calculate how the effective 
mixing angles @ij vary with the parameters of Hb and 5c p- Succintly put, we have 

Uf = Uf ({%} , {Ouj} , {Ami} ' {^iii ' ^cp, h, (t>b2, h^) = Uo ({G,,} , 6f) . (19) 

Up to this point, Uf is dependent on 14 free parameters. However, the standard mixing 
parameters Am^i, Am^i, 612, ^13 and 623 have been fixed by neutrino oscillation experiments. 
Additionally, in order to simplify the analysis, in the present work we have set all of the phases 
to zero, i.e. 5cp = = (i)2 = <i>3 = 0- Furthermore, also to simplify our analysis, we have made 
the eigenvalues of Hm proportional to those of Hb, at an energy of = 1 PeV, that is. 

Am?- ^ ^ 

bij = , (20) 

with A the proportionality constant. The upper bounds on the bij, Eq. ()16p . are satisfied for 
A < 10^. Thus we are left with only four free parameters to vary: A, 9bi2, Obis and 0623- 



3. Looking for extreme effects in the flavour ratios 

In the preceding Section, we saw that, in order for the energy-independent contribution to flavour 
transitions to be visible, we would need to use the expected high-energy astrophysical neutrino 
flux. As mentioned in Section [H the sources of this flux, e.g. active galaxies, are located at 
distances of tens to hundreds of Mpc, so that the appropriate flavour transition probability to 
be used is the average one, Eq. ([5]). 

If, at the sources, neutrinos of different flavours are produced in the ratios (j)^^ : (fP^ : 4>^, then, 
because of flavour transitions during propagation, the ratios at detection will be 

<t^a= (PM^ (21) 

/3=e,/^,T 

for a = e, /i, r. Evidently, the initial flavour ratios depend on the astrophysics at source, which 
is currently not known with high certainty. The standard initial flux [TO] considers that the 
charged pions created in proton-proton and proton-photon collisions decay into neutrinos and 
muons, which decay into neutrinos, too: 

7r+ /i+i/^ e+i/eI7^f^ , vr" e^Vei^pp^ ■ (22) 

Such process yields (p^ : (p^ : cp^ = 1 : 2 : 0. In the standard oscillation scenario, i.e. in the 
absence of the energy-independent contribution, plugging this initial flux into Eq. (I2ip . and using 
the best-fit values of the mixing angles, Eq. (fT2]) . yields equal detected fluxes of each flavour, 
i.e. (p'J'^ : (pf^ : (/.f « 1 : 1 : 1 

In a different production process [HI [l2l [13], the muons produced by pion decay lose all 
their energy before decaying, so that a pure-f^ flux is generated at the source, i.e. cfP^, : 
(p^ : (j)^ = : 1 : 0. In the standard scenario, this initial flux yields a detected flux of 
(pl^'^ : (p^^"^ : (p'^^'^ 0.22 : 0.39 : 0.39. Alternatively, a scenario of pure-z^e production 
(0g : : = 1 : : 0) through beta decay has been considered, e.g. in [11]. In this scenario, 
high-energy nuclei emmited by the source have sufficient energy for photodisintegration to occur, 
but not enough to reach the threshold for pion photoproduction. The neutrons created in the 
process generate through beta decay. For this initial flux, the resulting detected fluxes, in 
the standard oscillation scenario, are (pl^'^ : <p^J;'^ : cp^^'^ 0.57 : 0.215 : 0.215. In what follows, we 
will consider the possibility of observing the energy-independent contribution of Hb assuming 
these three initial flavour fluxes. 



Table 1. Standard values (no energy- independent contribution) of the detected fluxes 
[a = e, /i, r) and of the ratios R, S, for the three different initial fluxes considered in the 
text. The detected fluxes were calculated using the average flavour-transition probability in 
Eq. ([5]) with the central values of the mixing angles: sin2(0i2) = 0.304, sin2(0i3) = 0.01, 
sin^ (023) = 0.50. 
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We have defined the ratios 

^ = ^ ' ^=^' (23) 

and studied the effects of the new physics on them. The standard values of R and S, along 
with the standard values of the detected fluxes (pa, for each production scenario, are shown in 
Table [TJ Note that 5" = 1 for every production mechanism because the value of 6*23 used in the 
table was its best-flt value vr/4, which ensures yields equal detected fluxes of and Vr due to 
maximal mixing. 




1.0 1.5 2.0 2.5 3.0 3.5 4.0 

R 



Figure 1. Allowed regions of values of R and S assuming an initial flux of cp^ : (p^ : (p^ = 
1:2:0. The light blue region corresponds to the scenario when A = 100 (thus making the 
energy- independent contribution the dominant one) and 6bi2, Gbis, ^623 being allowed to vary, 
independently, within the range [0, vr] . The darker-coloured region, shown for comparison, is the 
allowed region of values obtained when A = (that is, energy-independent contributions turned 
off) and each of the standard mixing angles 9ij is allowed to vary within its la bound. A few 
notable points, marked A-D, are reviewed in Table [2j 



Figure [T] shows, in light blue, the allowed region of values of R and S when A = 100 and 
the mixing angles ^612, 6bi3 and 9b23 are allowed to vary, independently, within the range [0,7r], 
with the standard mixing angles set to their best-fit values, Eq. (I12p . This value of A already 
corresponds to the scenario where Hb dominates over Hm in the full Hamiltonian; higher values 
of A will not modify the shape of the allowed R — S region. For comparison, the darker-coloured 
region corresponds to the scenario where there is no energy-independent contribution (i.e. A = 0) 
and the standard mixing angles are allowed to vary within their la bounds. We have marked 
four notable points in the plot, A-D. Table [2] shows for each of them the values of the effective 
mixing angles Qij and the detected fluxes 

If we now consider the two other possible initial fluxes, (j)'^ : (j)^ : cj)^ = : 1 : and 1:0:0, set 
A = 100 to make Hi, the dominant contribution and allow each of the 9bij to vary within [0, tt], 
we obtain the R — S regions in Figure O The light purple and light brown regions correspond, 
respectively, to the initial fluxes 0:1:0 and 1:0:0 when A = 100. The light blue region is the 
same region that was shown in Figure [H assuming an initial flux of 1 : 2 : and A = 100. The 
dark blue, dark purple and dark brown regions correpond, respectively, to standard oscillations 
(A = 0), allowing the standard mixing angles 9ij to vary within their la bounds. 

Table 2. Notable points in Figure [TJ values of the effective mixing angles @ij and of the 
detected fluxes. 
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Because we are considering neutrinos that travel distances of tens or hundreds of Mpc, 
neutrino decay is a possibility. Assuming that the neutrinos decay into products that are not 
detectable (i.e. "invisible daughters" such as sterile neutrinos), then, following [T4j, the flux of 
flavour a at Earth will be 

'^"= E E'^°l[^o]/3/l[f^o]„jV^/^- E <^°l[f^oUl'l[f^o]./ , (24) 

fi=e,fi,T i /3=e,/i,T i(stable) 

where Tj is the lifetime of the i-th mass eigenstate in the laboratory frame. As explained in [14], 
this expression corresponds to the case where the decay has been completed when the neutrinos 
arrive at Earth. 

In a normal hierarchy, vi is the only stable state and so 

4>a = \[UoU\' E ^p\lUoW , (25) 

f3=e,fi,T 

while in an inverted hierarchy 1^3 is the stable state and 

'P» = \[UoW E ^'\Po]f3f , (26) 

f3=e,fi,T 

In Figure El decay to ui has been coloured red and decay to u^, green. We see that the R — S 
region for decay into vi does not intersect the region accessible with A = 100, assuming an initial 



flux of 1 : 2 : or : 1 : 0. It is, however, totally contained within the allowed region for 1:0:0, 
and thus there is a possibility of not being able to disentangle the effects of neutrino decay (in a 
normal hierarchy) from an energy-independent contribution, under this production model. The 
region allowed for decay into 1^3, on the other hand, is superposed to that of : 1 : 0, for values 
of R > 14 and 0.75 < S < 1.05. 
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Figure 2. Allowed regions of values of R and S for the different initial fluxes (jrl : (jr^ : 4w = 1 : 
2 : (blue), 0:1:0 (purple) and 1:0:0 (brown). Light blue, light purple and light brown 
correspond to A = 100, a situation where the energy-independent contribution is dominant 
in the flavour transitions. These regions were generated by varying each Ouj-, independently, 
within [0,7r], and fixing the standard mixing angles 9ij to their best fit values. Dark blue, dark 
purple and dark brown correspond to A = 0, that is, standard oscillations, without any extra 
contribution. These regions were generated by varying the Oij within their Icr bounds. Neutrino 
decay into a single lightest mass eigenstate has been considered, considering both a normal 
hierarchy (red) and an inverted one (green). These regions were also generated by varying the 
Oij within their la bounds. 



A notable feature of this plot is that it reveals that, while there are many opportunities of 
establishing the presence of an energy-independent contribution by measuring R and 5, it is 
a more difficult task to deduce from this measurement what the initial fiux was. An extreme 
example occurs in the region around R = 2, S = 0.2, where the three assumptions of initial 
flux result in the same prediction for R and S. In other words, if we measured values for R and 
S around this region, we could conclude (provided there are no other new physics effects) that 
there is an energy-independent contribution present, but we could not conclude what the initial 
flux was. On the other hand, if we measured a value of 4 < i? < 12, then (again, assuming 
there are no unaccounted effects), we could conclude both that there is an energy- independent 
contribution and that the initial flux was : 1 : 0. The same would occur if we measured i? < 1, 
S > 1.35: in this case, in addition to concluding that an energy-independent contribution is 
present, we could also conclude that the initial flux was 1:0:0. 

As evidenced in Figure [3l if we now make the contributions from standard, mass-driven, 
oscillations and energy- independent new physics comparable, effectively setting A = 1, the 
regions of allowed values of R and S do not change signiflcantly with respect to the case 



with A = 100, corresponding to dominance of the new physics. Hence, it is not necessary 
for the new physics effects to be dominant in order for them to affect the astrophysical flavour 
neutrino flux in a measurable way. In the A = 1 case, however, we recover some capability 
to disentangle neutrino decays to vi from the allowed region corresponding to the production 
mechanism 1:0:0. Additionally, for values of S < 0.1 and i? ~ 2, it seems to be possible to 
identify, in principle, the production mechanism to be 1 : 2 : 0. 
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Figure 3. Allowed regions of values of R and S for the different initial fluxes (f)^ : (jP^ : (ff^ = 1 : 
2 : (blue), 0:1:0 (purple) and 1:0:0 (brown). Light blue, light purple and light brown 
correspond to A = 1, a situation where the energy- independent contribution is comparable 
to the standard terms from mass-driven oscillations. Dark blue, dark purple and dark brown 
correspond to A = 0, that is, standard oscillations, without any extra contribution. 



Because S depends on the number of z^r's detected, and given that this is expected to be a 
low number (e.g. one every two years at IceCube |14]). the uncertainty on 5, using the present 
generation of neutrino telescopes, will be high. Figures [3] and [2] show estimated error bars at 
the points corresponding to standard oscillations using the best-fit values of the mixing angles, 
considering a 30% uncertainty on 5, due to the low statistics of Vr-, and a lower value of 15% 
for i?, dominated by the systematic error introduced by the atmospheric flux. Nevertheless, 
knowledge of only R could be enough to establish whether energy-independent new physics exists 
at the UHE scale or even to deduce what the initial was. For instance, we have seen that, if a 
value of R between 4 and 12 is measured then we could establish that there is energy- independent 
new physics and that the initial flux is : 1 : 0. 

4. Summary and conclusions 

We have shown that introducing an energy-independent contribution to neutrino flavour 
oscillations in the form of an additional term in the oscillation Hamiltonian affects the mixing 
angles, which, in turn, modify the flavour-transition probabilities. This contribution is motivated 
by scenarios of new physics such as flavour-dependent coupling to the gravitational field and 
violation of CPT symmetry. Given that the standard oscillation Hamiltonian has a 1/E 
dependence on energy, we have considered in our analysis the expected ultra-high-energy 
astrophysical neutrino flux, so that it becomes possible for the dominant term in the full 



oscillation Hamiltonian to be the energy-independent term. This flux, however, propagates 
over distances of tens or hundres of Mpc, so that we need to use the average flavour-transition 
probability and, more importantly, it makes sense to consider the possibility that neutrino decays 
occur. The modifications to the mixing angles result in a change in the expected flux of each 
flavour that arrives at Earth, which depends both on the initial flux of each flavour and on the 
oscillation probabilities. 

For the initial flavour fluxes, we have considered three possibilities: either neutrinos are 
produced by the decay of pions generated in proton-proton and proton-photon collisions, and 
of the muons produced by pion decay, which yields an initial flux of cp^ : (j)^ : (j)^ = 1 : 2 : 0; 
or, in a situation where the muons lose most of their energy before decaying, they are produced 
predominantly as direct products of the decay of pions, which yields a pure-f^ flux of : 1 : 0; 
or else they are produced by the beta decay of neutrons liberated in the photodisintegration of 
nuclei, yielding a pure-i^e flux of 1 : : 0. 

We have defined the ratios R = (j)f_i/(j)e and S = </>t-/</>^, and shown that large deviations 
from their standard values are possible, in all three neutrino production scenarios, given our 
current lack of knowledge about the intervening new physics parameters. Both when the effects 
of the new physics are dominant as well as when they are comparable to the standard oscillation 
terms, our ability to deduce, from a single measurement of R and S, the existence of an energy- 
independent contribution and the value of the initial flux is dependent on the particular values 
of R, S measured. Given that the number of i^'t-'s to be detected at current neutrino telescopes 
is very low (about one every two years at IceCube), however, the uncertainty on S will be quite 
high. In spite of this, knowledge of only R could be enough to establish the presence of the 
extra contribution or to deduce what the initial flux was, with particular sensitivity to a : 1 : 
initial flux. 
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